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A Hubbard model at half-filling on an anisotropic triangular lattice has been proposed as the 
minimal model to describe conducting layers of k — {BEDT — TTF)2X organic materials. The 
model interpolates between the square lattice and decoupled chains. The k — {BEDT — TTF)2X 
materials present many similarities with cuprates, such as the presence of unconventional metallic 
properties and the close proximity of superconducting and antiferromagnetic phases. As in the 
cuprates, spin fluctuations are expected to play a crucial role in the onset of superconductivity. We 
perform a weak-coupling renormalization-group analysis to show that a superconducting instability 
occurs. Frustration in the antiferromagnetic couplings, which arises from the underlying geometrical 
arrangement of the lattice, breaks the perfect nesting of the square lattice at half-flUing. The spin- 
wave instability is suppressed and a superconducting instability predominates. For the isotropic 
triangular lattice, there are again signs of long-range magnetic order, in agreement with studies at 
strong-coupling. 



I. INTRODUCTION 

The family of k — {BEDT — TTF)2X layered organic crystals exhibit many fascinating electronic properties0~i. 
There are similarities with the high- Tc cupratesQ. Competition between antiferromagnetic and superconducting in- 
stabihties, seen in the cuprates, also appears in the n — {BEDT — TTF)2X compounds. |-| 

The Hubbard model at half-filling on an anisotropic triangular lattice was proposed by McKenzieQ as a model 
of the conducting layers of k — {BEDT — TTF)2X. It is a simplification of a model first introduced by Kino and 
Fukuyamalj. Two hopping matrix elements are considered, ti between sites on a square lattice and ^2 between 
next-nearest-neighbors along one of the two diagonal directions (see Fig. |]). 




FIG. 1. Anisotropic triangular lattice with two hopping amplitudes ti and t2. 



The case ti = <2 thus corresponds to the isotropic triangular lattice. The model interpolates between the square 
lattice (^2 = 0), which has been the subject of many studies in the context of high-Tc cuprate superconductors and 
which has a spin-wave instability at half-filling, and completely decoupled chains {ti = 0) for which we have an exact 
Bethe ansatz solution by Lieb and Wucl. The Hamiltonian is given by: 



H = -ti^{cl''cja + h.C.)-t2 (c^^Cj^T +/I.C.) C/o^(71i| - - ^) 
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where < ij > denotes nearest- neighbor pairs of sites on the square lattice and << ij >> are next-nearest- neighbor 
pairs along one of the two diagonal directions as shown in Fig. [|. _ _ 

Values for the hopping amplitudes, obtained from several quantum chemistry calculationsuij and for different 
anions X, are in the range ti > t2', that is, somewhere between the square and the isotropic triangular lattices. The 
non-interacting part of the Hamiltonian can be easily diagonalized and the corresponding Fermi surfaces at half-filling 
are depicted in Fig. |[ 
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FIG. 2. Fermi surface of non-interacting electrons for different values of the hopping amphtudes. The number on the top 



of each graph is the value of the quantity t^jitY + i2), ranging from (square lattice) to 1 (decoupled chains), 
potential n is varied to ensure that the system is at half-filling. 



The chemical 



In the next section we present rewilts from a weak-coupling renormalization-group analysis on this model. We use 
the approach of Zanchi and SchulalJ. The Fermi surface is divided into 16 patches. One-loop RG flow equations for 
all the possible scattering processes involving the patches on the Fermi surface are numerically integrated towards 
the low-energy limit. All subleading non-logarithmic contributions, arising from six-point functions generated during 
the mode elimination, are included. We conclude with a discussion of our results and make a comparison with results 
obtained in the strong-coupling limit. 



II. RESULTS 



The Hubbard model on the pure square lattice (ti = 0) and with repulsive on-site interaction C/q > has been 
extensively studied. Weak-coupling renormaliaatipn-group analyses of the couplings, resolved into discrete patches 
along the Fermi surface, have been performedlJ^ll^. At half-filling, since the Fermi surface is fully nested (Fig. ||), an 
antiferromagnetic (AF) spin-wave instability develops. As the system is doped slightly away from half-filling, there 
is a crossover to a superconducting BCS regime and the dominant Cooper pairing symmetry is (1^2 _y2. As shown in 
Fig. Isl we have reproduced these resultslHnlj. 
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FIG. 3. Results for the square lattice (t2 = 0). At half-filling {fj, = 0), spin- wave instability occurs, but as the system is 
doped away from half-filling, the perfect nesting is lost and there is a crossover to a BCS instability. 



Introducing non-zero t2 offers a different way of favoring a BCS instability over AF. Perfect nesting is eliminated 
once the hopping t2 between next-nearest-neighbors sites along one of the two diagonal directions is turned on (Fig. 
|l|) . In contrast to the square lattice where both AF spin- wave and BCS couplings increase during the renormalization- 
group transformations with V^^ diverging faster than y-^cs^ sufficiently large t2 there is a crossover to a regime 
where the BCS processes eventually dominate, signaling a superconducting instability. Furthermore, because the 
Fermi surface is imperfectly nested, the growth in both types of couplings weakens. Further increasing t2 eventually 
destroys the nesting of the Fermi surface altogether and both types of divergences are suppressed. Three cases 
illustrating these crossovers are shown in Fig. 0. 



2 



40.0 




0.0 2.0 4.0 6.0 8.0 

scaling parameter I 



FIG. 4. At half-filling, spin-wave instability occurs for t2 — 0, but as t2 increases, the BCS instability wins over. Finally, as 
t2 is increased further, the nesting is destroyed and both divergences are suppressed. 

Going further to the isotropic point, our weak-coupling RG analysis no longer shows BCS instabihties. In Fig. |^ 
the dominant AF and BCS channels are compared. Neither channeLshows strong divergences, but the AF channel is 
significantly larger than the BCS channel. Slave-boson calculationstJ find a Mott-Hubbard metal-insulator transitioa 
to occur at the relatively large value of Uc = 7.23t, but also find a metallic phase with incommensurate spiral orderll3 
at C/ < Uc- Our weak-coupling results seem to indicate the onset of re-entrant long-range antiferromagnetic order. 
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FIG. 5. Flow of dominant AF and BCS channels for the isotropic triangular lattice. 



AF ordering tendencies again disappear as t2 is increased further, beyond ti. This is as expected, since decoupled 
Hubbard chains do not exhibit AF order. 



III. CONCLUSION 



Hubbard models on the square lattice have been extensively studied in the context of high- Tc superconductivity. 
We have reproduced the well-known result that there is an AF instability at half-filling. Also doping the system away 
from half filling induces a crossover to a BCS regime with dx^-y2 pairing symmetry, as expected. More importantly, 
we have shown another way of triggering a BCS instability. Keeping the system at half-filling, but introducing the 
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diagonal hopping t2 (as shown in Fig. |T]), eUminates perfect nesting. Corresponding magnetic frustration kills the 
spin density wave, and Cooper pairing can dominate. This result suggests that superconductivity may occur in a 
model of strongly-correlated electrons that interpolates between the square and isotropic triangular lattices. 

At the isotropic point, there are signs of re-entrant antiferromagnetic long range order. In the large-f7o limit, at 
half-filling, the Hubbard model can of course be mapped to the spin-1/2 Hcisenberg antiferromagnet which is known 
to have long-range AF order on the triangular lattice. Furthermore, away from the isotropic point, our weak-coupling 
RG analysis continues to agree qualitatively with results obtained for the pure Heisenberg system with two exchange 
couplings Ji and J2. Thejshase diagram of the Heisenberg model has been studied via a straight 1/S expansionllj, 
series expansion methodstll, and a large-N treatmenttSl. All three methods find two regions of long-range order: near 
the limit of a square lattice (J2 = 0) and near the isotropic point (J2 = Ji). It is remarkable that our weak-coupling 
analysis agrees with these results. 
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